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■ ABSTRACT 

'-^ ' GRB 061007 is the brightest gamma-ray burst (GRB) to be detected by Swift and is 

i accompanied by an exceptionally luminous afterglow that had a l/-band magnitude 

Q ■ < 11.1 at 80 s after the prompt emission. From the start of the Swift observations the 

S-j \ afterglow decayed as a power law with a slope of ax = 1-66 ± 0.01 in the X-ray and 

c/3 . a op t = 1-64 ± 0.01 in the UV/optical, up to the point that it was no longer detected 

above background in the optical or X-ray bands. The brightness of this GRB and the 
, similarity in the decay rate of the X-ray, optical and 7-ray emission from 100 s after 

*h . the trigger distinguish this burst from others and present a challenge to the fireball 

' model. The lack of a cooling or jet break in the afterglow up to ~ 10 5 s constrains 

any model that can produce the large luminosity observed in GRB 061007, which 
we found to require either an excessively large kinetic energy or highly collimated 
outflow. Analysis of the multi-wavelength spectral and high-resolution temporal data 
taken with Swift suggest an early time jet-break to be a more plausible scenario. This 
must have occurred within 80 s of the prompt emission, which places an upper limit 
on the jet opening angle of 8j = 0.8°. Such a highly collimated outflow resolves the 
energy budget problem presented in a spherical emission model, reducing the isotropic 
equivalent energy of this burst to E^ orr = 10 50 ergs; consistent with other GRBs. 
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1 INTRODUCTION 

Gamma-ray bursts (GRBs) are by far the most energetic 
and luminous astrophysical phenomena in the Universe, re- 
leasing isotropic equivalent energies of the order of 10 48 — 
10 54 ergs in 7-rays alone on timescales of 10 -3 — 10 3 s. 
The favoured mechanism to produce such a large amount 
of energy on such a short timescale is the creation of 
a black hole, either from a merger of two compact ob- 
jects or the gravitational core collapse of a massive star. 
There appear to be two distinct populations of GRBs, dis- 
tinguished by the duration of their prompt emission and 



their spectral properties. The most robust model for long, 
soft GRBs, which typically have prompt emiss ions lasting 
> 5 s, is the collapsar model (|Wooslevl IT993). in which 
a GRB results from the energy released during the death 
of a massive star. Evidence supporting this model includes 
the large number of lon g GRBs located in regions of high 
star-f ormation rate (e.g. Tanvir et al.l |2004| ; iFruchter et al.l 
2006), and the association of several near by GRBs with 
super novae, e.g. GRB 980425 /S N1998bw jKulkarni et al] 
ll99Sfi . GRB 030329 /S N2003dh (iHiorth et all 120031 ) and 
GRB 060218/SN2006aj (|Campana et alj|2006h . 
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There is evidence to suggest that collimated outflows re- 
duce the amount of energy emitted by a GRB from isotropic 
emission by ~ 2 orders of magnitude, where the beam- 
corrected energy is clustered at 10 50 — 10 51 ergs (Frail et 
al. 2001; Bloom, Frail & Kulkarni 2003). Bipolar jets are 
a natural consequence of accretion onto a black hole and 
are observed in many high-energy astrophysical phenomena 
(e.g. AGN, pulsars, micro-quasars), where the ejecta take 
the easiest exit path, along the rotational axis. 

A signature of highly relativistic jets is an achromatic 
break in the broadband afterglow light curves, caused by 
the expansion of the radiation beaming angle as the bulk 
Lorentz factor, F, decreases. Once the Lorentz factor has 
decreased such that 6j = 1/F, where 9j is the jet opening 
angle, the rate of decay increases, resulting in a break in the 
light curve. 

Swift is a rapid-response mission with multi-wavelength 
capabilities designe d specifically to stu dy the early time 
evolution of GRBs |Gehrels et alj|2004 ). The spacecraft is 
equippe d with three telescope s; the Burst Alert Telescope 
(BAT; iBarthelmv et al1l2005h . which detects 7-rays from 
~ 10 up to several hu ndred keV, the X-Ray Telescope (XRT; 
iBurrows et al.l 120051 ') which covers the 0.3-10 keV energy 
range, and the Ultr a- Violet and Optical Telescope (UVOT; 
iRoming et all f2005). which can observe in six colours and a 
broadband white band filter, covering the wavelength range 
from 1600 A to 6000 A. 

Swift's slewing capabilities provide GRB afterglow light 
curves from ~ 60 s after the trigger, and well-sampled X-ray 
and optical light curves have revealed far fewer achromatic 
breaks than had been expected prior to launch (Ghirlanda, 
Ghisellini & Lazzati 2004; O'Brian et al. 2006). The jet open- 
ing angle, 9j, increases with tj as 8j tx t^ 8 , and thus the 
lack of observed jet breaks may be the result of large average 
jet opening angles. 

GRB 061007 is an example of a burst that showed no 
evidence of a break in the X-ray or UV/optical afterglow 
from 80 s after the prompt emission, until it was no longer 
detected above background, at ~ 10 6 s. Prompt observa- 
tions by ground based observatories measured the redshift 
of this burst to be z=1.26 (Osip, Chen & Prochaska 2006; 
Jakobsson et al. 2006). At this redshift, GRB 061007 had 
one of the highest isotropic equivalent energies ever seen, 
releasing ~ 10 54 ergs in 7-radiation alone, w ith an optical 
after glow comparable to that of GRB 990123 (|Akerlof et al.l 
1 19991 ), with V < 11.1 mag at T + 71.7 s. 

In this paper we describe in detail the 7-ray, X-ray and 
optical observations taken by Swift. From the high resolu- 
tion early time data we constrain the type of environment in 
which the burst occurred and determine the conditions that 
distinguish this burst from others. In section [5] we describe 
the data reduction and analysis techniques used followed by 
a discussion on the physical parameters and emission mech- 
anisms involved that could reproduce our observations in 
section 13.113.2.31 We also analysis the properties of the cir- 
cumburst medium of GRB 061007 in section 13.41 using the 
afterglow spectral energy distribution (SED), and a sum- 
mary of our results is given in section [4] Throughout the 
paper all errors quoted are la unless otherwise specified 
and the temporal decay rates, a, and spectral indices, 0, 
are defined such that Ft tx t~ a and F v cx respectively. 
Luminosities and the energy release are calculated assuming 



0.3 
0.2 
0.1 









6 


-p 






0) 






X) 





4 








CO 





2 


w 






-p 






un 


























6 


CD 






■+-> 






ra 





4 








d 





2 























E-i 






< 






CQ 







0.2 
0.1 




^1 1^^!^.:^^^^ 




Figure 

the four 



-50 50 100 150 

Time since BAT trigger (s) 

1. Prompt emission light curve of GRB 061007 shown in 
BAT energy bands. 



a standard cosmology model with Ho = 70 km s 1 Mpc 1 , 
D,m = 0.3 and fi A = 0.7. 



2 OBSERVATIONS AND ANALYSIS 
2.1 BAT 

The BAT on board Swift triggered on GRB 061007 on 7 th 
October, 2006 at 10:08:08 UT, resulting in an immediate 
slew to point the narrow field instruments in the direction 
of the burst. From here on we shall refer to the time of 
the trigger as T. The refined BAT position for GRB 061007 
is RA = OS^OS'm.S 3 , Dec = -50°29'47.4" (J2000) with 
an uncertaint y (systematic and stat istical) of 0.9' at 90% 
containment ( Markwardt et akll2006h . 

The time-interval over which 90% of the 15-150 keV 7- 
radiation is emitted is T90 = 75 ± 5 s. The BAT light curve 
consists of three large and distinct peaks (see Fig. [TJ with 
weak emission detected out to ~ 150 s after the trigger. 
The time lag between the peak emission in the 50-100 keV 
and the 15-25 keV energy bands was 16lg ms, and between 
the 100-350 keV and 25-50 keV bands the lag was 21±g ms 
(private communication with J. Norris). 
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2.2 XRT 

The first XRT observation was a centroiding image taken 
80.4 s after the BAT trigger, and a bright uncata- 
iogued source was found at RA = 03 h 05 m 19.5 s , Dec = 
-50°30'01.9" wit h a 90% containment radius of 3.5" 
(|Vetere et aljbood ). 

At 86.6 s after the BAT trigger the XRT began tak- 
ing data in Window Timing mode, in which data are stored 
with a 1.7 ms time resolution and 1-dimensional imaging. 
At T + 2000 s the count rate fell to below ~ 1 count s" 1 and 
the automated sequence switched to Photon Counting (PC) 
mode, in which the time resolution is 2.5 ms and full imaging 
and spectroscopic resolution are available. The X-ray after- 
glow was detected for 10 8 s after the BAT trigger, before it 
fell below the detection threshold. During the first 2200 s of 
observations the count rate of the burst was high enough to 
cause pile-up in both WT and PC mode data. To account 
for this effect, the WT data were extracted in a rectangular 
40x20- pixel region with a 9 x 20-pixel region excluded from 
its centre. The size of the exclusion region was determine d 
following the procedure illustrated in lRomano et al.l(|2006al ). 
The standard correct ion for pile-up was ap plied to the PC 
data, as described bv lVaughan et all ||2006T ). Due to a large 
amount of pile-up in the first 30 s of WT observations, these 
data were not used in our spectral analysis. 

The 0.3-10 keV light curve can be represented well by 
a power law with decay index a = 1.66 ± 0.01 from T + 80 s 
until at least T + 10 6 s, when it reached the detectability 
threshold with XRT (see Fig. [5} • We also determined the 
earliest times that a cooling break or jet break could fit the 
data with a goodness of fit within 3<r of the best-fit power 
law model. To do this we applied a broken power law fit 
modelled on the cooling break, for which cY2 = ot\ + 0.25, 
and on a jet break, where 02 = Qi + 0.75. In the former 
model t Va > 10 5 s before the change in \ 2 from a power 
law fit was more than 3cr level, and in the latter model this 
corresponded to at tj > 2.65 x 10 5 s. 

A power law fit to the WT spectrum, where the 
Galactic c olumn density is fixed at Nh = 2.13 x 
10 cm" " 2 (|Dickev fc Lockmanlll990l) . gives a best- fit host 
galaxy column density at a redshift of z — 1.26 of (5.1 ± 
0.3) x 10 21 cm -2 for a spectral index /3 X = 0.99 ± 0.02. 



2.3 UVOT 

The UVOT began observing GRB 061007 71.7 s after the 
BAT trigger, at which point it took a 9 s settling exposure 
in the U-filter. This was followed by a 100 s exposure in the 
white light filter (A ~ 1600 - 6500 A) and a 400 s exposure 
in the U-band. After this the automated sequence rotated 
twice through the UVOT filters taking a series of short ex- 
posures (V,UVM2,UVW1,U, B, white, UVW2) (10 s for 
B and white filters and 20 s for the rest). The central wave- 
lengths to the optical and UV filters are 5460 A, 4370 A 
3460 A, 2600 A, 2260 A and 2030 A for the V, B, U, UVW1, 
UVM2 and UVW2, respectively. A further 400 s U-band 
exposure and 100 s white filter were taken, followed by an- 
other series of rapid rotations around the filter wheel up to 
2260 s after the burst. All observations up to this time were 
taken in event mode, which has a time resolution of 11 ms, 
after which time image mode exposures were taken. The re- 



fined position of the afterglow is RA = 03 h 05 m 19.6 s , Dec 
= -50°30'02.4" to a certainty of 0.5" (J2000). The Galac- 
tic extinction along this line of sight is E(B — V) = 0.021 
(Schlegel, Finkbeiner & Davis 1998). 

For all short exposures (20 s) and exposures taken later 
than T + 2260 s photometric measurements were obtained 
from the UVOT data using a circular source extraction re- 
gion with a 6" radius for the V , B and U optical filters and a 
12" radius for the three UVOT ultra-violet filters to remain 
compatible with the current effective area and zeropoint cal- 
ibrations 1 . The background was measured in a source-free 
region near to the target using an extraction radius of 20". 

The first U-band settling exposure and first white-band 
exposure were affected by coincidence loss to a degree that is 
outside the photometrically calibrated range of UVOT; for 
a U-band settling exposure, complete saturation within a 
standard 6" aperture implies that this afterglow is brighter 
than V — 11.1 for the whole of the 9 s exposure. During the 
first white band exposure the source is so bright that rela- 
tive photometry can be obtained using a very large aperture, 
which includes the wings of the point spread function, where 
coincidence loss is not a problem. The exposure was divided 
into 5 s bins, and source counts were obtained from an ex- 
traction region of radius ~ 26". To calibrate this photom- 
etry, observations of the white dwarf standards GD50 and 
HZ2 were taken with UVOT and the white filter an October 
18 th 2006. From these observations, we determined the rela- 
tionship between observed count rate within the 26" aper- 
ture and corrected count rate within a standard 6" UVOT 
source aperture, which we assume to be linear. To apply this 
calibration to GRB 061007 we apply constant offset to this 
relation to account for the different background count rate 
between the GRB 061007 observations and the white dwarf 
standard count rate observations. This offset was chosen to 
match the mean white-band count rate in the saturated ex- 
posure to the white filter light curve extrapolated from later 
times; we caution that calibration of this first white exposure 
is crude and thus the apparent steepening of the light curve 
towards the end of the exposure may not be real. However, 
the overall behaviour of a steadily decreasing count rate over 
the course of the exposure is not in doubt. For the first 400 s 
U-band exposure a higher resolution light curve was created 
by splitting it into 4 s bins, and the second into 20 s bins, 
consistent with the rest of the exposures during that time. 

We created a single UV/optical light curve from all the 
UVOT filters in order to get the best measurement of the op- 
tical temporal decay. To do this the light curve in each filter 
was individually fitted to find the corresponding normalisa- 
tion, and this was then used to renormalised each light curve 
to the U-band. A fit to the combined UV and optical light 
curve, from T + 85 s up to T + 2.5 x 10 4 s, is best fit with a 
power law decay with temporal index a opt = 1.64±0.01. The 
combined light curve shows no apparent colour evolution, as 
can be seen in Fig. [2] 

The combined UVOT light curve is shown next to the 
0.3-10 keV and 15-350 keV light curve in Fig. H The mag- 
nitude along the left axis corresponds to the U-band, and 
the count rate along the right axis applies to the X-ray and 
7-ray light curves. 



1 http: / /heasarc . gsfc . nasa.gov / docs/heasarc /caldb / swift /docs/ uvot / 
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Figure 2. GRB afterglow light curve shown in all six UVOT lenticular filters, in the X-ray (0.3-10 keV) and in the 7-ray (15—350 keV). 
All light curves have been shifted vertically with respect to the V band filter for clarity. The B, U, UVW1, UVM2, UVW2 and white 
filter have all been normalised to the V band light curve. The left hand axis corresponds to the V-band light curve, and the count rate 
axis on the right applies to the X-ray and 7-ray light curves. The best fit decay indices to the combined UVOT light curve, the X-ray 
light curve and the BAT light curve are ajjvOT = 1-64 ± 0.01, cexRT = 1-66 ± 0.01 and ct^AT = 1-54 ± 0.12, respectively. 



2.4 Observations By Other Facilities 



The prompt emission from GRB 061007 also triggered 
Konus-Wind and Suzaku, and the afterglow was monitored 
by ground based telescopes, which carried out both spec- 
troscopic and photometri c observations. GRB 061007 was 
meas ured to be at z=1.26 l|Qsip et al. 2006; Jako bsson et al.l 
2006), for which the rest-frame isotropic energy release and 
maximum luminosity in the 20 keV-10 M eV is Ej so ~ l.Ox 
10f* erg, and L iso ~ 1.8 x 10 53 erg s" 1 ijGolenetskii et all 
l2006h . The optical afterglow was also detected by ROTSE 
at 26.4 s after the BAT tr igger at a magnitude of 13.6 mag 
in an unfiltered exposure (jRvkoff fc RuioDakarnl l2006). and 
the Faulkes Telescope South measured an fl -band magni- 
tude of 7? = 10.15T0.34 142 s after the trigger |Bersier et all 
2006), consistent with the UVOT observations and indica- 
tive of a significant brightening during the first two minutes. 
The Swift observations from T + 80 s show that by this time 
the afterglow had entered a constant power law decay (see 
Fig. [2}. Radio observations were also performed 1-1.24 days 
after the prompt emission, which resulted in no detection 
of a radi o source with a flux measu rement of 1 ± 36 [iJy at 
19 GHz |van der Horst fc RollfioOr^ l. 



2.5 Multi- Wavelength Light Curve Analysis 

At ~ T + 100 s, after the last flare in the BAT light curve 
(Fig.[T|, the BAT light curve appears to decay smoothly and 
monotonically. A power law fit to the BAT light curve from 
~ T + 100 s provides an acceptable fit (x 2 = 32 for 29 dof ) 
with a best-fit decay index of olbat = — 1.54±0.12, which is 
consistent with the decay observed in both the UV/optical 
and X-ray light curves. Furthermore, a power law spectral fit 
to this last part of the BAT light curve, beyond ~ T + 100 s, 
gives a best fit energy spectral index /3bat(> T + 100 s) = 
—0.80 ± 0.09, in fairly good agreement with the spectral 
indices determined for the X-ray and optical afterglows. This 
supports a scenario in which the 7-emission at this epoch is 
generated by the same radiation mechanism as the X-ray 
and UV/optical afterglow. 



2.6 The Spectral Energy Distribution 

Using the full range of data available with Swift a SED was 
produced using optical, X-ray and 7-ray data, which were 
normalised to an instantaneous epoch at 600 s after the BAT 
trigger, at which time there are data available from all three 
instruments, and good colour information for UVOT. For 
each of the UVOT lenticular filters the tool uvot2pha (vl.l) 
was used to produce spectral files compatible with xspec 
(v12.2.1), and the count rate in each band was normalised 
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Figure 3. The combined 7-ray, X-ray and UV/optical spectral energy distribution of GRB 061007 at T + 600 s (see text) with best-fit 
models for each corresponding dust extinction curve shown; SMC (dashed), LMC (solid) and Galactic (dotted). Open circles are data 
points at wavelength A < 1215 A in the rest frame, and therefore not used in the spectral fitting. The rest-frame positions of the Galactic 
absorption feature at 2175 A the beginning of the Lyman-a forest at 1215 A and the Lyman-a break at 912.5 A are indicated. 



to the estimated value at T + 600 s determined from the 
best-fit models to the light curves. The XRT spectrum was 
produced from WT mode data in the 0.5-10 keV energy 
range, omitting the first 30 s of data, which suffered from 
pile-up. This was scaled to the corresponding count rate at 
T + 600 s. For the BAT spectrum, data were taken from 
the time interval T + 100 s - T + 500 s, during which the 
BAT light curve decayed at a rate consistent with X-ray and 
UV/optical afterglows and had a similar spectrum. Beyond 
500 s the 7-ray emission became indistinguishable from the 
background. This spectrum was normalised to T + 600 s in 
the same way as the XRT spectrum. The SED is shown in 
Fig.E 

The SED was fit with a power law emission spectrum 
and two dust and gas components to model the Galactic and 
host galaxy photoelectric absorption and dust extinction. 
The column density and reddening in the first absorption 
system were fixed at the Galactic values (see sections 12.21 
and I2.3|l . The second photoelectric absorption system was 
set to the redshift of the GRB, and the neutral hydrogen 
column density in the host galaxy was determined assum- 
ing solar abundances. The dependence of dust extinction on 
wavelength in the GRB host galaxy was modelled on three 



extinction laws, taken from observations of the Milky Way 
(MW), the Large Magellanic Cloud (LMC) and in the Small 
Magellanic Cloud (SMC). The greatest differences observed 
in these extinction laws are the amount of far UV (FUV) ex- 
tinction, which is greatest in the SMC and least in the MW, 
and the strength of the 2175 A absorption feature, which is 
most prominent in the MW and negligible in the SMC. To 
model these extinction laws the parameterisations from IPeil 
(1992) were used, where R v = A V /E(B - V) = 3.08,2.93, 
and 3.16 for the Galactic, SMC and LMC extinction laws 
respectively. 

At a redshift of z — 1.26 the beginning of the Lyman- 
a forest is redshifted to an observer-frame wavelength of 
~ 2580 A, which falls squarely on the UVW1 filter, the 
reddest of the UV filters. The unknown equivalent width 
of Lyman-a absorption from the host galaxy of the GRB 
and the variation of the Ly-a forest along different lines 
of sight makes it difficult to model this additional source 
of absorption in the SED. We, therefore, exclude the three 
data points affected by Lyman-a from our spectral analysis. 
The results from our spectral analysis are summarised in 
Table[T] The fitted value of the spectral index depends on the 
amount of optical/UV extinction, which in turn is affected 
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Table 1. Power law spectral fit results to the spectral energy distribution of GRB 06f007 



Spectral Model 


Host galaxy 
Nh 1 (10 21 cm" 2 ) 


Rest-Frame 

visual extinction, Ay 


Po,X 


x 2 

(dof) 


Null- hypothesis 
probability 


SMC 
LMC 
MW 


3 66 +0-29 

r on+0.33 

O,au -0.32 

7 62+ - 81 
'■ DZ -0.77 


0.39±0.0f 
0.66 ±0.02 
1.13 ±0.10 


0.90 ± 0.005 
0.98 ± 0.007 

-1 -,0+0.04 
1 ' lo -0.03 


398 (308) 
347 (308) 
1042 (308) 


4 X 10~ 4 

0.063 

0.000 



Equivalent hydrogen column density assuming solar abundances at the GRB host galaxy. 



by the dust extinction model used. Similarly, the amount of 
soft X-ray absorption fitted to the model will affect the best- 
fit spectral index. Therefore, although the X-ray absorption 
and optical/UV extinction are independent components in 
the fit, they are coupled because they both depend on the 
spectral index. 

A fit to the optical, X-ray and 7-ray data, with an LMC 
extinction law for the GRB host galaxy, gives an acceptable 
fit, with x 2 ~ 347 for 308 degrees of freedom (dof), im- 
plying a null-hypothesis probability of p = 0.063. Models 
with an SMC and MW extinction law at the host galaxy 
were rejected with 99.96% confidence and 100% confidence, 
respectively. 

3 DISCUSSION 

GRB 061007 is one of the brightest bursts to be ob- 
served by Swift, with an isotropic energ y released in 7- 
rays that is second only to GRB 050904 (Tagliafcrr i et all 
120051 ) and an optical flux compa rable to that observed in 
GRB 990123 (lAkerlof et alj|l999h . The similarity of the af- 
terglow of GRB 061007 in the flux decay of 7-ray, X-ray and 
UV /optical bands, and the lack of features in X-ray and op- 
tical light curves make this GRB highly unusual. Optical 
afterglows typically decay as power laws, although usually 
with decay indices ~ 1.1 (Zeh, Klose & Kann 2006). Indeed, 
prior to the launch of Swift GRB 061007 could have been 
categorised a dark burst. In most GRB X-ray afterglows, 
the light curve is characterised by several discrete segments, 
where the canonical shape is made up of three, and some- 
times four, power law segments with dif ferent decay indices 
l|Nousek et al.ll2006l ; Izhang et alj|2006ah . 

The lack of a break in the light curves provides a con- 
straint on the frequency of the characteristic synchrotron 
frequency, v m and the peak flux density, _F^ imax . From the 
response to the white filter it can be seen that around half 
the counts correspond to the V and B filter, and the lack 
of a break in the first white filter exposure, therefore, indi- 
cates that at T + 80 s v m was no bluer than the B-band 
filter. This corresponds to an upper limit on the character- 
istic synchrotron frequency v m < 7.5 x 10 14 Hz. From the 
UV/optical data we also determine that at T±80 s the flux 
density in the U-band, corresponding to v — 5.5 x 10 14 Hz, 
is F v = 730 mjy, after correcting for absorption at the host 
galaxy. 

With these constraints in place we investigate two pos- 
sibilities that could explain the temporal and spectral prop- 
erties of GRB 061007, where the fundamental difference is 
in the effect of collimation on our observations. In the first 
model we assume that a jet break does not occur before 



T + 2.65 x 10 5 s, based on the constraints we obtained in 
section 12.21 and the observed emission is, therefore, indis- 
tinguishable from isotropic emission. We then look into a 
model in which the jet break occurs before our first X-ray 
and UV/optical observations at T + 80 s such that assump- 
tion of isotropic emission no longer applies. 

3.1 Spherical Expansion Model 

The relation between a and /3 for synchrotron emission at 
different intervals in the synchrotron spectrum are defined 
by a set of closure relations, both for a constant density 
mediu m (Sari, Piran fc Nara yan 1998) and a wind density 
profile (|Chevalier fc Lilll999l) . Using these we find our ob- 
servations to be in closest agreement with a constant den- 
sity circumburst medium, where the UV/optical and X-ray 
bands are in the same slow-cooling regime (v m < v < v c ). 
In this case, the temporal and spectral flux decay indices 
are related as a = 3/3/2. For a value of (3 X = 0.99 ± 0.02, 
we expect a — 1.49 ± 0.03, which is shallower than the ob- 
served decay rate of a = 1.65±0.01. However, this deviation 
from the relation is small, and could be accommodated if the 
density profile is slightly deviated from an absolute constant 
medium. In this regime /3 = (p — l)/2, where p is the elec- 
tron energy index. Thus, for /3x = 0.99 ± 0.02 the electron 
energy index is in the range 2.94 < p < 3.02. We therefore 
take p = 3.0. These relations hold as long as the character- 
istic synchrotron frequency, v m , lies below the optical band, 
i.e. v m < 7.5 x 10 14 Hz, at T±80 s, as discussed previously. 
There is also the additional condition that the cooling fre- 
quency, u c , must lie above the peak X-ray frequency up to 
at least T + 10 5 s, before which a cooling break would be 
detected in our X-ray observations (see section 12. 2p . This 
corresponds to the requirement that v c > 6 x 10 17 Hz at 
T+10 5 s. 

For a constant density medium in which synchrotron 
radiation is the dominating emission mechanism, the char- 
acteristic synchrotron frequency, the cooling frequency and 
the peak density flux are given by the following equations 
from lZhang et all (|2006bh 

00^ m 1 5tj/'P~2\ 2 , 1 ,1/2 1/2 2t-i1/2 ,—3/2 /-t\ 

u m = 3.3x10 Hz^-— j-J (l + z)' e B e cE k ; 52 t d (1) 

v c = 6.3X10 12 Hz(l+z)- 1/2 (l+Y)- 2 e B 3/2 E i ;^ 2 n- 1 t d " 1/2 (2) 

F^max = 16 mjy(l + z)D2 8 2 eB /2 E ki5 2n 1/2 , (3) 

where es is the fraction of energy in the circumburst mag- 
netic field, e e is the fraction of energy in electrons, Ek, 52 is 
the kinetic energy in the GRB in units of 10 52 , td is the time 
since the prompt emission in units of days, n is the density 
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of the surrounding circumburst medium, Y is the Compton 
parameter, and D28 is the luminosity distance in units of 
10 28 cm. 

For a spectral index (3 = 1, vF v is constant for u m < 
v < u c . At T + 80 s F y = 730 mjy at v = 5.5 x 10 14 Hz, 
and therefore i/ m F„, max = 730 x (5.5 x 10 14 ). By using this 
together with Eq. [T] and [3] we obtain: 

3.4 x 10~ 6 

71 = l T e l E? (4) 

In Eq. [5] we apply the condition that at 10 5 s after the 
prompt emission v c must still be greater than 6 x 10 17 Hz, 
and substituting in Eq. [4] gives 



0- + Y) 



-2 1/2 4 j-,5/2 



fpK ci > 0.52 



(5) 



In the case of the characteristic synchrotron frequency 
we know that v m < 7.5 x 10 14 Hz, and this reduces Eq. [T]to 



4 /2 £e<52 < 1-7 X 10- 



(6) 

This, and the inequality given in Eq. [5] provide an upper 
limit to the value of eb in terms of I?f 52 , such that 



e B < 3.1 x 10- 19 (l + y)- 4 ^|, 52 



(7) 



The above equation bounds us to extreme values of E k ,52 
and e B . For E k ^ 2 = 10 3 , e B < 3.1 x 10^ 10 (1 + V) -4 , which 
is an unrealistically low value of eg. Even a value as small 
as es = 10~ 6 requires a huge amount of kinetic energy to 
meet the condition that E k , 52 > 1.5 x 10 4 (1 + Y) 4/3 . The 
amount of kinetic energy assumed is an isotropic equivalent 
value which would, thus, be smaller if it were corrected for 



collimation. A 3<r lower limit of ti 



2.65 x 10 5 s can be 



imposed as the earliest time that a jet break could have 
occurred, tj is the time of the jet break measured from the 
onset of the afterglow, which we take to be equal to T. A 
kinetic energy E kfi2 = 1.5 x 10 4 (1 + F) 4/3 and e B = 10~ 6 
would require e e = 0.012 and n = 49 cm -3 for the conditions 
given above to be satisfied. From this we estimate the jet 
opening angle using 



= 0.2 



3 / 8 / n 
V j 



1/8 



(8) 



.1 + zJ V£fc,52 

l|Panaitescu fc Kumar! 120001 ). where tj t a is in units of days. 
For the physical parameters assumed, a jet opening angle 
for GRB 061007 would need to be 6j > 6.3° to produce a 
jet break later than T + 2.65 x 10 5 s. This would provide 
a beam-corrected kinetic energy of E k orr > 9 x 10 53 ergs, 
which is three orders of magn itude greater than usual bursts 
(IFreedman fc Waxmanl feOQl). The energetics required by 
this model are unreasonably large and, therefore, an unlikely 
solution to explain the properties of GRB 061007. 



3.2 Highly Collimated Models 

The large amount of energy required to power a model with 
isotropic emission leads us to investigate a scenario in which 
the emission from the GRB is highly collimated, where we 
still assume the circumburst environment to be of constant 
density. A very narrow jet opening angle could reduce the 
energy required to explain our observations by several orders 
of magnitude, where tj must be prior to our first afterglow 
observations in order to account for the lack of a break in 



the X-ray and UV/optical light curves at later times. If we 
assume the afterglow begins at the time of trigger, this cor- 
responds to a time tj < 80 s. Within the premise that the 
outflow is highly collimated we consider two cases in which 
the cooling frequency is either above or below the observ- 
ing frequencies. In these two regimes p and /3 are related 
by /3 = (p — l)/2 and f3 = p/2, respectively, corresponding 
to p ~ 3.0 or p ~ 2.0. After the jet break the flux decays 
as F t cx t~ v (Sari, Piran & Halpern 1999), where p is ei- 
ther 2 or 3, which is steeper than the observed decay rate 
of a ~ 1.65. However, this does not account for continual 
energy injection, either from the central engine itself, or late 
time refreshed shocks that slow down the decay of the af- 
terglow. The shallow decay phase observed in many XRT 
light curves is interpreted to be the result of energy injec- 
tion, and although this phase typically lasts for 10 4 s, there 
are also examples of GRBs with evidence of energy injection 
for prolonged periods of time, and with very smooth decay 
slopes (e.g. Romano et al. 2006; Huang, Cheng & Gao 2006). 
Taking this into account, an explanation whereby continual 
energy injection maintains the afterglow and slows down its 
decay over a long period of time is not unreasonable. 



3.2.1 u m < v < v c 

For a laterally expanding jet with constant density circum- 
burst medium and v m < v < u c , the relation between the 
temporal decay of the GRB afterglow, a, and the spectral 
index, (3, post jet-break is given by 



a=(l + 2/3)--(l 



9)08 + 2), 



(9) 



|Panaitescu et all 120061) where q is t he energy injection 
param eter such that L = Lo(t/to)~ q , l)Zhang fc Meszarosl 
l200lf ). For our observed temporal and spectral indices, q = 
0.3 and p = 3, as previously. Prior to the jet break Eq. I ll'il 
still apply, but in order take into account the effect of energy 
injection, the kinetic energy requires an additional factor of 
{t/tj) ■ The time of jet break and the jet opening angle 
are related as 9j cx t 3//8 . To keep 9j as large as our observa- 
tions will allow we take tj = 80 s. At T + 80s the additional 
factor of (t/tj)^ 1- ^ is unity and there is no change in Eq. [T} 
El Thus at T + 80 s, Eq. 0] and Eq. [6] stand. For Eq. [2] the 
constraint that v c > 6 x 10 17 at T + 10 5 s means that the 
evolution of v c post jet-break has to be taken into account. 
If v Cj is the frequency of v c at T + 80 s, 



a nc 1 n!9 !/ 2 4 TP 5 / 2 

v Cj = 4.05 x 10 e s ' e e E k ' 



(10) 



However, we are concerned with where v c lies at T + 10 5 s. 
After a jet break v c cx 7 -4 t -2 , where 7 cx t~( 2 + 9 )/ 6 is the 
Lorent z factor of the electron s in the circumburst environ- 
ment (|Panaitescu et al1 l2006). The cooling frequency after 
T + 80 s is given by 



1(9-1) 



(11) 



where jj is the electron Lorentz factor at T + 80 s. Using 
the condition that v c > 6 x 10 17 Hz at T+ 10 5 s, for q = 0.3 
Eq. [H] reduces to 



(i + y)-4 /2 4< 5 2 2 >o.4i 



(12) 
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In the same way that we determined the inequality given 
in Eq. [7] we combine the above equation with Eq. [6] and 
determine that 



es < 5.0 x 1CT 19 (1 + Y)- 4 E 3 k , 



(13) 



For £ fe ,52 = 5x 10 4 , e B < 6.2 x 1(T 5 (1 + V)" 4 , and if we take 
this upper limit as the value of £s, e e = 3.1 x 10~ 3 (1+Y) 1 ' 2 . 
Although the assumed kinetic energy is large, this corre- 
sponds to the isotropic equivalent value and the beam- 
corrected kinetic energy is a few orders of magnitude smaller 
than this. For p — 3, the inverse Compton parameter Y is 
proportional to [j?(e e /es)] 1/ ' 2 , where ?j=min[l, (v-m/vc) 1 ], 
and for these values of e e and es, Y < 0.3. For values, of 
£k,52 = 5 x 10 4 , eg = 6.2 x 10~ 5 , and e e = 3.1 x 10~ 3 
the particle density in the circumburst environment is n = 
0.08 cm -3 , which, using Eq. [8] gives an opening angle of 
6j = 2.0 x 10~ 3 rad. This corresponds to 6j = 0.11°, which 
gives a beam-corrected kinetic energy E%° rr = 10 51 and 7- 
ray energy E^° rr = 2 x 10 48 . The jet opening angle that 
we estimate assumes the h ydrodynamics as described by 
iPanaitescu fe Kumarl l|2000h . and variations on the hydro- 
dynamics, such as that discussed in I Sari et alj (Il999l ). result 
in a different coefficient in Eq. [8] However, these differences 
are of order unity, and do not affect the validity of the model 
discussed. 

This model is, therefore, able to satisfy our observations 
without the requirement of an excessively large kinetic en- 
ergy. However, the opening angle i s a factor of at lea st 30 
narrower than previously observed dBloom et al.l [20031 ') and 
the inferred 7-ray energy is extremely low for a long GRB, 
prompting us to consider other options. 



3.2.2 v > max[i/ m , v c ] 

Keeping in place the condition that tj < 80 s, we investigate 
a model in which the X-ray and UV/optical bands are above 
v m and v c . In this model our observational constraints on 
the values of v m , v c and F„, max only apply at T + 80 s. 
Since we assume a jet break to occur at tj — 80 s, the 
factor of (i/i,)' 1-9 ' in the kinetic energy is unity, and q is 
therefore not of importance in our calculations. For v c above 
the observing band, /3 = p/2, and therefore p = 2, in which 
case the dependence on the minimum Lorentz factor of the 
shock accelerated electrons 7 m changes from 



_ (P-2) mp 
(p - 1 ) m c 



to 



7m 



In 



(-) 



m e 



(14) 



(15) 



where e m and sm are the minimum and maximum energy 
of the shock accelerated electro ns, and V is the bulk Lorentz 
factor in the shocked medium l|Pe'er fc Zhand[2006l ). v m oc 
7m , and therefore the expression for v m given in Eq. [1] has 
to be modified such that [(p — 2)/{p — l)] 2 is replaced by 
[ln(e m /eM)] 2 - The value of \ii(em /em) is not well determined 
due to the unknown magnetic field in the upstream of the 
GRB shock, although it generally ranges between 5 and 10 
|Li fc Waxmanl 120061 ). We therefore take \n(e M /e m ) = 7, 
such that Eq. [T] becomes 



We assume the same jet break time as before (tj = 80 s) 
whereby Eq . [2l and 1161 reduce to 



z~b ' Kl n(l + Y)-' < 5.45 
and 

4 2 *lE\%< 2.1x10- 
In the regime where v > max[i/ m , v c ], 
F — F i/ 1|/2 i/ 



,l/2„(p-l)/2 I/ -p/2 



(17) 

(18) 
(19) 



By taking the measured value of F v = 730 at v = 5.5 x 
10 14 Hz, and using Eq. [3] [T7] and 1181 the expression for F„ 
reduces to 



e e £fc, 6 2(l + rr 



3.65 



(20) 



a -7 iniS tt / 1 1 \l/2 1/2 2^1/2 ,-3/2 

v m = 6.7 x 10 Hz(l + z) e B ' e c E k ; 52 t d 



(16) 



From these equalities we then determine that a kinetic en- 
ergy of £ fei52 = 10 3 would require e e = 3.7 x 10~ 3 (1 + Y), 
e B < 0.25(1 + F)~ 4 andn > 0.05(l + y) 4 . For a particle den- 
sity of n — 10 4 cm -3 in the GRB local en vironment, which 
is a t ypical density for a molecular cloud (|Reichart fc Price] 
2002), the jet opening angle would be 8j — 0.8°, a factor 
of a few smaller th an that determined f or some other GRBs 
(e.g. GRB 980519; iNicastro et"aHll999h . This gives a beam- 
corrected kinetic energy of = 9.7 x 10 50 . 

This model provides an alternative scenario in which a 
highly collimated jet is the cause for the large flux observed 
in GRB 061007, but where a larger jet opening angle and 
more standard beam-corrected 7-ray can satisfy our obser- 
vations, with respect to the previous model. 



3.3 Synchrotron Self-Compton Emission 

For completeness we look into a scenario in which syn- 
chrotron self-Compton (SSC) emission contributes to the 
emission. In the case where the X-ray emission is domi- 
nated by SSC radiation, and the UV/optical emission is pro- 
duced by synchrotron radiation, both observed bands would 
need to be in the fast-cooling regime, with v opt > v c > v m 
and vx > Vc C > v m \ m order for the decay slope and 
spectral index in the X-ray and UV/optical emission to 
be consistent. In this case the spectral index is given by 
P = p/2, which would give a value of p = 2.0. However, 
the rate of synchrotron optical afterglow is expected to be 
— (2 — 3p)/4, which is inconsistent with our observations of 
a= 1.65 ±0.01. 



3.4 The Circumburst Medium 

From our spectral analysis we determine the dust extinction 
in the circumburst medium to be most consistent with a 
wavelength dependence similar to that of the LMC. This is 
consistent with previous studies of GRB environments (e.g. 
Fruchter et al. 2006, Kann, Klose & Zeh 2006; Schady et al. 
2006), which also find GRB hosts to be irregular galaxies, 
and is also in agreement with the collapsar model, which re- 
quires a sub-solar metallicity progenitor star, with an upper 
limit of Z Q < 0.3 (e.g. Hirschi, Meynet & Maeder 2005; 
Woosley & Heger 2006). 

The best-fit parameters for the X-ray column density 
and dust extinction local to the GRB are Nh = (5.3 ±0.3) x 
10 21 cm _2 and Av = 0.66 ±0.02, and this gives a gas-to-dust 
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ratio of N H /A V = (8.0 ±0.32) x 10 21 . However, N H was de- 
termined assuming solar abundances, which are unlikely to 
be correct for an irregular, Magellanic-type galaxy. Correct- 
ing the X-ray column measurement to a metallicity of 1/3 
solar, appropriate for an ISM which is similar to that of the 
LMC, we find that the gas-to-dust ratio is a factor of four 
larger than that measured in the SMC or LMC. The gas-to- 
dust ratio measured in the local environment of GRB 061007 
is also larger than the avera ge of those deter mined for a set 
of previous Swift GRBs (jSchadv et alj|2006h by a factor of 
~ 1.5. This could be the result of dust destruction in the 
environment of the GRB. 

The lack of evolution observed in the soft X-ray absorp- 
tion later than ~ 80 s after the prompt emission implies that 
the burst must have completely ionised the gas in its im- 
mediate environment within this time scale. Given the large 
column density that is measured from the SED we, therefore, 
infer that the absorbing medium we observe is far enough 
from the burst that it was not completely ionised in the first 
few 100 s of seconds after the initial explosion. The probabil- 
ity of there being an absorption system in the line of sight 
from an unrelated, intervening galaxy that is of sufficient 
optical depth to contribute a measurable amount of soft X- 
ray absorption is ver y small (|Q'Flahertv fc Jakobsenlll997l ; 
ICampana et al.l [20061 ). leading us to conclude that the ab- 
sorbing material measured in the spectru m of GRB 061007 
is loca l to the burst. From the results in iPerna fc Lazzatil 
(|2002h we estimate the absorbing medium to lie at approx- 
imately ~ 10 pc from the GRB. The large value of Nh is 
comparable to that of giant molecular clouds, from which 
we infer that GRB 061007 was embedded in a dense, star- 
forming region, or giant molecular cloud, in common with 
many other GRBs (e.g . iDai fc Lul Il999l ; (Wang et al.1 12000| ; 
iReichart fc Pricell2002h . 

3.5 Comparing Between Models 

It is interesting to note that with E lt i so = 10 54 ergs in the 
10 keV-20 MeV energy range, and a rest-frame peak en- 
ergy of E p = 902 keV, GRB 061007 satisfies th e Amati 
relation, given by E p tx E^ so (I Amati et al.ll2002h . It also 
satisfies a correlation between the peak isotropic luminos- 
ity, L Pt i SO , and spectral lag in the prompt emission, Ti ag , 
whereby L p> i sa oc Tj™ -1 . Both these correlations are sat- 
isfied by a large fr action of long GRBs (e.g. lAmatil [20061 ; 
iGehrels et aT1l2006l ), suggesting that GRB 061007 has the 
same class of progenitor as other long GRBs. 

The consistency observed in the spectral and tempo- 
ral afterglow behaviour in the X-ray and UV/optical bands, 
and in the late time 7-ray band, implies that the emission 
observed comes from the same spectral segment. If we take 
the afterglow to be in the slow cooling regime and assume 
a jet break to occur at tj > 2.65 x 10 5 s a large amount of 
kinetic energy, of at least E^ orr = 9 x 10 53 is required. The 
implausibly large amount of energy required by this model 
makes it an unsatisfactory scenario to explain the properties 
of GRB 061007. 

We also investigate a second scenario in which a highly 
collimated outflow produces a jet break before T + 80 s. 
This requires a jet angle no larger than 0.8°, which satisfies 
the observational requirements with a lower energy budget 
and more reasonable set of parameters. The upper limit on 



the jet opening angle of 0.8° is a factor of 4 smaller than 
that previously determined in GRBs (Berger, Kulkarni & 
Frail 2003), although still physically possible. Furthermore, 
the need to observe the jet break in order to measure 8j 
introduce selection effects, and the smaller and larger end 
of the current jet opening angle distribution will be under- 
estimated. For the very smallest of jet opening angles to be 
measured, very early time data of the afterglow is needed, 
preferentially covering several energy bands. 

A consequence of a very early jet break time is the faster 
evolution of V into a non-relativistic phase, which will pro- 
duce a break in the light curve as the decay rate slows down. 
However, F cx i _1//2 , and therefore the GRB outflow will 
not enter a non-relativistic phase until t nr = tj0~ 2 . For 
8j < 0.8°, this corresponds to t nr > 4 x 10 5 s, which is not 
in conflict with our observations. 

For a jet break earlier than T + 80 s our observations 
are valid for a spectral regime in which v m < v < u c , and 
v c < v, where in the latter a larger jet opening angle is 
obtained by assuming a density n = 10 4 cm -3 , consistent 
with that observed in molecular clouds. The dependence of 
the circumstellar particle density on es, e e and Et,52 in the 
former of these models does not allow for such high densities, 
for which we estimate n = 0.08 cm -3 . Our spectral analysis 
on the afterglow resulted in a large column density, which 
supports a scenario in which the GRB is embedded within 
a dense molecular cloud or active star-forming region. This 
therefore favours a model in which v c < v, where the number 
density is inferred to be larger. 

In the case of a uniform jet, where the energy and 
Lorentz factor are constant across the jet, the afterglow after 
a jet break should decay as t~ p . In both the post jet-break 
models that we consider the predicted decay slope is steeper 
than the observed value of a = 1.65±0.01. A possible expla- 
nation is continual energy injection, which would maintain 
the energy in the afterglow and slow down the rate of de- 
cay. The requirement for steady energy injection over the 
entire light curve is of some concern, but seems plausible 
on the basis of many other examples of energy injection in 
Swift afterglow observations. Evidence of prolonged energy 
injection in the form of plateaus in X-ray light curves is 
observed in a large fraction of Swift GRBs, although typi- 
cally up to 10 3 -10 4 s after the prompt emission. In the case 
of GRB 061007 a constant amount of injection is needed 
for at least an order of magnitude longer than this, during 
which no break is observed in the light curve. This is not 
to say that the central engine needs to be active for this 
length of time, but simply that it is active for the duration 
of the prompt emission and ejects shells with a large distri- 
bution in t he Lorentz factor. A range in T from a few tens 
to ten (e.g. iRees fc Meszarosl [l99sl ; ISari fc Meszarosl l200Ch 
can provide a source of energy to the afterglow for a long 
time after the central engine has stopped through ongoing 
refreshed shocks. Nevertheless, an energy injection mecha- 
nism that remains extremely constant for such a long time 
requires very specific conditions in which the energy injec- 
tion mechanism dominates the afterglow during the entire 
period of observations. 

Alternatively, the slow decay rate of the afterglow in the 
context of post jet-break could also be the result of a jet- 
edge effect in which there is no sideways expansion. In such 
a case the afterglow decay index would steepen by 3/4 when 
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r = 1 /0j , and would conti nue to decay at this rate due to the 
lack of lateral expansion ijMeszaros fe Reedll999h . In such 
a scenario there is no need for energy injection, although 
additional conditions would be required to prevent the jet 
from expanding further, such as a high magnetic field within 
the jet. 

A possibly more natural form of non-laterally expanding 
jet is a structured jet or outflow with an energy distribution 
that is a decreasing function of angular distance from the jet 
axis, 80, and a narrow, high-energy core at the axis. A jet 
break is, then, observed when the edge of the uniform inner 
core becomes visible to the observer. In order for a jet break 
to occur before T+80 s the viewing angle, therefore, needs to 
be close to the jet axis, which is not that unreasonable given 
the brightness of GRB 061007. At the time of jet break the 
smoother distribution in energy, e, compared to a uniform 
jet results in a smaller change in the decay slope that will 
be shallower than the 3/4 increase in decay index observed 
in a uniform jet break. For a power law distribution of E, 
where E cx 8~ k , the smaller the value of k the shallower 
the afterglow light curve will be after the jet-break. In the 
case of v m < v \, s < v c , where p — 3, a decay rate after a 
jet break with index a = 1.65 will occur in structured jet 
with an E distri bution index fc = 0.3, and for v c < v \, s we 
require k = 1.4 (|Panaitescull2005l ). 

In the context of post-jet break emission, a decay slope 
of a = 1.65 can, therefore, be accounted for by a non- 
laterally expanding jet as well as continual energy injection 
through refreshed shocks. We, therefore, find an early jet 
break before T + 80 s can satisfy the conditions imposed 
by our observations without requiring a large energy budget 
and with a reasonable set of physical parameters. 



4 CONCLUSIONS 

In this paper we have presented multi-wavelength spectral 
and temporal data of GRB 061007 taken with Swift; an in- 
credibly bright GRB with unusually smooth panchromatic 
decay. Using these data we explored the circumstances that 
could produce such a bright burst, and so determined the 
conditions that distinguish GRB 061007 from other, more 
standard GRBs. On many levels we find GRB 061007 to 
be typical of others, suggesting that the same progenitor 
model for most long GRBs is also applicable to GRB 061007. 
The surrounding dust has properties most consistent with 
an LMC-like environment, and the circumburst material is 
likely to be dense, suggesting that it is embedded in a large 
molecular cloud or a region of active star formation. Fur- 
thermore, its consistency with both the Amati and timelag- 
luminosity relations, which apply to a large fraction of long 
GRBs, indicate that the source of energy in GRB 061007 is 
no different than that of most other long GRBs. 

The enormous brightness of GRB 061007 is, therefore, 
unlikely to be due to a different progenitor model, but the 
result of extreme values for a few defining characteristics. 
Two inherently different models can satisfy our observations 
of GRB 061007, where the main difference lies in the time of 
the jet break; either before or after our observations. Both 
require certain extreme conditions to account for the bright- 
ness of this burst, either in the form of a large energy budget, 
or in the collimation of ejecta into a very narrow jet open- 



ing angle. However, given the brightness of GRB 061007, it 
should not be a surprise that some of defining properties are 
non-standard. The large amount of kinetic energy required 
in a late-time jet break model lead us to favour a model in 
which the GRB 061007 is produced by the concentration of 
energy into a highly narrow jet < 0.8°; a condition that can 
more easily be met than the production of large amounts of 
energy. 

The brightness of GRB 061007 across the electromag- 
netic spectrum makes it ideal to explore the conditions 
surrounding the production of GRBs, and its uniqueness 
provides further insight into the range in properties that 
GRBs have. The small opening angle that we determine for 
GRB 061007 has important implications on the true range 
in jet opening angles, and this could provide a clue as to 
the cause for the large fraction of Swift GRBs with no clear 
signs of jet breaks. 
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